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ABSTRACT 


The  NPS  computer  simulation  model  was  modified  to  study 
sputtering  by  molecular  ions.   The  simulations  were 
performed  on  a  Cu(lll)  surface  for  normally  incident  O2 
ions  at  0.5,  1.25,  2.5,  and  5.0  kev/ion.   The  molecule's 
angular  orientation  was  both  set  at  specific  values  and 
randomly  determined.   The  normalized  energy  distributions 
of  sputtered  atoms  for  two  energy  ranges,  0-20  and  0-3 
ev/sputtered  atom,  and  the  sputtering  yield  ratios  were 
investigated.   The  0-20  ev/sputtered  atom  energy 
distributions  show  distinct  peaks  at  1.8  ev  which  are 
insensitive  to  incident  energy,  energy  density,  and 
molecular  orientation.   Simulations  indicate  that  there  may 
be  fine  structure  in  the  0-3  ev/sputtered  atom  range  with  a 
cascade  interaction  unique  to  molecular  sputtering.   For  a 
number  of  molecular  orientations  which  produce  very 
different  cascade  overlaps,  the  sputtering  yield  ratios  do 
not  correlate  with  the  degree  of  overlap.   The  energy 
density  of  a  cascade  is  not  an  accurate  predictor  of  the 
sputtering  yield  ratio. 
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I.   INTRODUCTION 

A.   OVERVIEW 

Sputtering,  the  ejection  of  atoms  from  a  solid  by 
energetic  ion  bombardment,  is  caused  by  the  collision  of 
incident  ions  with  atoms  in  the  surface  and  near-surface 
layers . 

Sputtering  was  first  observed  by  Grove  [1]  in  1852. 
While  conducting  experiments  with  electrical  discharges  in 
a  low-pressure  gas,  he  noticed  that  cathode  material  was 
deposited  on  the  inner  wall  of  the  discharge  tube.   Thus, 
sputtering  gained  its  reputation  as  a  contaminating 
side-effect  of  electrical  discharge  processes.   This 
reputation  persisted  for  nearly  100  years  and  scientific 
interest  in  the  sputtering  phenomenon  was  largely  directed 
toward  its  elimination  rather  than  its  understanding. 
However,  major  scientific  achievements  in  the  20th  century, 
such  as  high  voltage,  high  power  radar  transmitter  tubes, 
focused  attention  on  sputtering.   Thereafter,  sputtering 
research  followed  two  parallel  courses:   First,  the 
development  of  a  basic  theoretical  understanding  of  the 
physics  of  sputtering,  and  the  acquisition  of  experimental 
data.   Second,  exploitation  of  this  knowledge  for  practical 
purposes . 


Sputtering  as  a  bulk,  macroscopic  process  was 
"understood"  and  early  applied  to  a  limited  degree.   For 
example,  the  development  of  pumps  based  on  ion  trapping 
processes  and  the  production  of  atomically  clean  surfaces. 
However,  very  little  true  understanding  of  the  physics  of 
sputtering  was  achieved  through  early  experimentation  or 
theoretical  models. 

B.   EARLY  EXPERIMENTATION 

Experimental  results  proliferated  but,  in  general,  the 
data  were  plagued  by  non-reproducibility  caused  by 
inconsistent  experimental  conditions  [2,3].   Penning  and 
Moubis  [4]  first  demonstrated  conclusively  the  effect  of 
ambient  pressure  on  sputtering  experiments.   They  also 
coupled  a  magnetic  field  to  their  discharge  apparatus  to 
reduce  back-diffusion  of  sputtered  material  to  the  surface 
[4].   Their  experiments  produced  the  first  consistently 
reproducible  results  for  energies  greater  than  500  ev. 

Surface  cleanliness,  the  elimination  of  oxide  layers  or 
adsorbed  gases,  is  a  critical  factor  in  sputtering 
experiments.   Arifov,  et .  al.,  [5]  demonstrated  in 
controlled  experiments  the  dependence  of  the  sputtering 
yield  on  surface  contamination.   Yonts  and  Harrison  [6] 
presented  evidence  indicating  that  surface  recontamination 
from  background  gases  is  a  significant  factor  in 


quantitative  sputtering  yield  measurements  and  developed 
criteria  for  "clean"  surface  experiments. 

C.   DEVELOPMENT  OF  ANALYTIC  THEORY 

Stark  [7,8]  first  extensively  discussed  the  concept  of 
an  individual  sputtering  event  on  an  atomic  scale  and 
developed  two  competing  models:   He  originated  the  hot-spot 
model  which  proposed  the  ejection  of  surface  material  from 
a  microscopically  small  region  due  to  evaporation  at  high 
local  temperatures  generated  by  ion  impact  [7].   His 
second,  the  collision  model,  described  the  sputtering  event 
as  a  series  of  binary  collisions  initiated  by  a  single 
incident  ion.   He  applied  the  conservation  laws  for  elastic 
collisions  to  this  model  and  qualitatively  described  the 
energy  dependence  of  the  yield  curve  [8].   Stark  obtained 
limited  agreement  with  experiments  at  energies  less  than  1 
kev. 

Kingdon  and  Langmuir  [9]  further  developed  the  binary 
collision  momentum  and  energy  model  to  describe  the 
observed  desorption  of  a  thorium  monolayer  from  a  tungsten 
surface  by  ion  bombardment. 

In  1926,  von  Hippel  [10]  formally  presented  the 
hot-spot  theory  of  sputtering  which  was  further  developed 
by  Townes  [11]  in  1944.   Townes '  theory  correlated  well 
with  some  contemporary  experimental  results,  but  the 
agreement  was  subsequently  shown  to  be  fortuitous. 


Keywell  [12]  developed  a  theory  based  upon  the 
moderation  of  neutrons,  which  may  be  taken  as  a  precursor 
of  the  collision  cascade  concept  to  follow.   Harrison  [13] 
developed  the  first  direct  statistical  theory  of  sputtering 
which  bridged  the  gap  between  the  Kingdon  and  Langmuir 
momentum  transfer  theory  and  Townes'  hot-spot  theory.   A 
significant  feature  of  the  Keywell  and  Harrison  theories  is 
the  application  of  probability  concepts  as  expressed  by 
collision  cross-sections. 

At  this  point  in  the  evolution  of  sputtering  research, 
the  vast  majority  of  experimental  data  were  based  upon 
sputtering  from  polycrystalline  surfaces  [3].   Sputtering 
was  assumed  to  be  isotropic  with  the  angular  distribution 
of  sputtered  material  following  a  cosine  distribution. 
Effects  due  to  the  underlying  crystal  structure  were  not 
considered. 

A  major  turning  point  in  sputtering  theory  was  the 
experimental  demonstration  of  lattice  effects  in 
monocrystal  sputtering.   Wehner  [14]  showed  that  sputtering 
from  single  crystals  proceeded  in  preferred  directions 
which  correlated  with  the  directions  of  closest  packing  in 
the  particular  crystal.   This  sputtering  anisotropy 
completely  contradicted  the  hot-spot  model  and  pointed  to 
momentum  transfer  as  the  key  ejection  mechanism.   Silsbee's 
focusing  collision  model  [15]  quickly  followed  as  a  partial 
explanation  for  preferential  ejection.   In  addition, 
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experimental  results  showed  the  sputtering  yield  to  be 
strongly  sensitive  to  the  crystallographic  direction  of  the 
incident  ion  beam.   The  subsequent  discovery  of  channeling 
[16-20],  and  the  elucidation  of  assisted  focusing  and 
replacement  collision  (previously  known  as  crowdion)  models 
[21]  showed  that  even  a  single  momentum  transfer  theory 
must  be  an  aggregate  of  many  different  mechanisms. 

D.   COLLISION  CASCADE  THEORY 

Sigmund  [22,23]  developed  an  extension  of  the 
statistical  theory  which  displaced  the  focusing  models.   It 
is  based  upon  a  linearized  form  of  the  general  Boltzmann 
transport  equation.   Sigmund' s  model  treats  an  amorphous 
solid  with  a  surface  binding  energy  potential  barrier.   An 
incident  ion  strikes  a  surface  atom  and  creates  a  collision 
cascade  via  binary  collsions.   An  atom  in  the  cascade 
sputters  if  it  is  headed  outward  from  the  surface  and 
overcomes  the  surface  potential  barrier.   In  doing  so,  the 
atom's  velocity  vector  is  refracted  away  from  the  surface 
normal  by  the  attractive  surface  binding  energy.   Signifi- 
cant agreement  with  high  energy  experimental  results,  that 
is,  incident  ion  energy  much  greater  than  the  surface 
binding  energy,  has  been  achieved.   However,  it  has  since 
been  shown  that  the  major  result  of  the  Sigmund  theory  fol- 
lows immediately  from  a  dimensional  analysis  of  the  systeml 
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Thompson  [24]  and  Robinson  [25]  have  taken  similar 
approaches  while  accounting  for  focusons  and  channeling, 
respectively,  and  have  obtained  comparable  results. 

The  primary  utility  of  the  Sigmund  and  Thompson 
theories  has  been  to  provide  experimentalists  with  simple 
formulae  to  apply  to  a  wide  variety  of  sputtering 
experiments.   However,  these  theories  fail  to  adequately 
predict  sputtering  yields  at  low  energies  [22,23,26].   The 
high-energy,  binary  collision  model  for  the  collision 
cascade  does  not  account  for  the  predominance  of  low-energy 
sputtered  atoms  from  near-surface  layers. 

In  general,  the  collision  cascade  can  only  be 
reasonably  understood  as  a  many-body,  or  multiple 
interaction  (MI),  process.   However,  the  MI  model  cannot  be 
solved  analytically. 

E.   COMPUTER  SIMULATION  MODELS 

Gibson,  Goland,  Milgram  and  Vineyard  (GGMV)  [27]  first 
demonstrated  the  power  and  feasibility  of  the  use  of 
computer  simulation  techniques  to  model  radiation  damage 
experiments.   Their  computer  program  builds  a  finite, 
three-dimensional  monocrystal  lattice  of  copper, 
accelerates  a  lattice  atom  in  a  specific  direction  with 
specific  energy,  and  solves  the  resultant  MI  cascade  by 
applying  numerical  methods  to  Newton's  equations  of  motion. 
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In  this  particular  case,  GGMV  studied  the  generation  of 
lattice  vacancies  and  interstitials . 

Robinson  and  Oen  [16,  17]  applied  GGMV's  prototype  to 
ion  range  studies  and  demonstrated  ion  channeling,  which 
was  subsequently  verified  experimentally  [18,  20]. 

Harrison  expanded  upon  GGMV's  concept  and  designed  a 
simulation  model  to  specifically  study  sputtering  [28-32] 
with  the  flexibility  to  simulate  a  wide  variety  of 
experimental  conditions.   Harrison's  early  results  showed 
that  neither  Focusons  nor  Sigmund's  or  Thompson's  models 
accounted  for  the  majority  of  sputtered  atoms  (see  also  the 
work  of  Lehmann  and  Sigmund  [33]).   Simulations  show  that 
the  vast  majority  of  sputtered  atoms  originate  in  the 
surface  layer,  and  the  results  agree  reasonably  well  with 
experimental  yield  data.   Simulations  augment  laboratory 
results  and  allow  the  isolation  and  manipulation  of 
variables  not  readily  accessible  in  the  lab  [34,  35]. 

F.   INTERATOMIC  POTENTIAL  FUNCTIONS 

A  primary  thread  of  uncertainty  that  runs  through  each 
sputtering  model  and  theory  is  the  form  of  the  interatomic 
potential  function  [3,  35,  36].  No  single  interatomic 
potential  has  yet  been  deduced  which  fits  every  sputtering 
sputtering  experiment.  Simulations  are  particularly  well- 
suited  for  isolating,  modifying,  and  observing  the 
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sensitivity  of  the  sputtering  event  to  the  nature  of  the 
interatomic  potential  [35,  37,  38]. 

G.   NON-LINEAR  CASCADE  EFFECTS 

Recent  experimental  results  indicated  that  the  analytic 
sputtering  theories  do  not  adequately  predict  the 
sputtering  yield  for  heavy  ions  bombarding  heavy  targets 
[39-41].   A  spike  theory  has  emerged  which  states  that 
linear  cascade  theory  breaks  down  (becomes  non-linear)  when 
the  average  energy  per  atom  in  the  cascade  approaches  the 
surface  binding  energy  of  the  target  material  [42].   The 
spike  model  emphasizes  random  ejection  effects  from  the 
region  of  high  energy  density.   Kelly  [43]  proposes  a 
thermal  sputtering  theory,  a  Renaissance  of  the  hot-spot 
theory,  to  explain  the  non-linear  effects  in  the  context  of 
the  spike. 

These  theories  reinforce  the  multiplicity  of  sputtering 
phenomena,  but  continue  to  fall  short  of  providing  a 
comprehensive  analytic  model. 

H.   MOLECULAR  SPUTTERING 

Bombardment  by  molecular  ions  is  well-suited  for 
investigating  non-linear  effects  since  the  yield  from  a 
single  ion  can  be  compared  to  the  yield  from  a  molecular 
ion  of  the  same  element  with  the  same  energy  per  ion.   Few 
published  results  of  molecular  sputtering  exist  [39-41,  44] 
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The  available  results  do  demonstrate  non-linear  effects  but 
do  not  confirm  the  abnormally  high  yield  enhancement 
addressed  by  the  spike  model  [42]. 

It  is  not  yet  clear  whether  the  individual  ion  cascades 
resulting  from  the  molecule's  dissociation  upon  impact 
overlap  to  significantly  increase  cascade  energy  density, 
or  whether  molecular  orientation  affects  the  sputtering 
mechanisms . 
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II.   OBJECTIVES 

Few  results  exist  for  molecular  sputtering.   Bader,  et . 
al,  [44]  conducted  low-energy  (<  8  kev)  experiments  with 
N+  and  N^  on  several  metals.   Their  primary- 
objectives  were  to  study  the  effects  of  the  upper 
atmosphere  on  satellites  and  other  space  vehicles  and  to 
investigate  erosion  rates  within  ion  propulsion  drive 
systems.   Recent  experiments  in  non-linear  cascade  effects 
use  high-energy  (100-200  kev/ion)  polyatomic  molecules  as  a 
convenient  method  of  increasing  primary  energy  density  in 
near-surface  layers  of  polycrystalline  targets  [39-42],   No 
other  experiments  have  been  conducted  to  study  the  effects 
of  molecular  sputtering,  per  se. 

The  primary  objective  of  this  thesis  is  to  study  the 
sputtering  of  a  copper  monocrystal  by  low-energy  (0.5  -  5.0 
kev/ion)  diatomic  molecules  using  a  computer  simulation 
model. 

More  specifically,  the  following  areas  will  be 
studied: 

1.  Comparison  of  energy  and  angular  distributions  of 
sputtered  atoms,  and  absolute  yields  for  molecular 

(0j)  and  ion  (0+)  bombardment. 

2.  Effects  of  molecular  orientation  upon  the 
sputtering  yield. 
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3.  Comparison  of  Oj/204"  yield  ratios  with 
published  results  for  Nj/2N+  [44]. 

4.  Comparison  of  0j/20+  and  ArJ/2Ar+ 

yield  ratios  for  indication  of  yield  enhancement  in 
accordance  with  non-linear  cascade  theory  [39-42]. 

5.  Effects  of  mass  ratio  (M^/M2)  and  ion  size 
(ion-atom  potential)  upon  yield  ratios. 

6.  Differences  in  incident  ion  scattering  for 
molecular  and  ion  bombardment. 

7.  Cascade  overlap. 
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III.   THE  MODEL 

A .   BACKGROUND 

QRAD,  a  full-lattice  simulation  computer  program,  is 
designed  to  model  the  dissipation  of  an  incident  ion's 
momentum,  using  classical  mechanics,  in  a  solid  monocrystal 
[28-32].   See  references  [34]  and  [35]  for  detailed 
discussions  of  the  theory  and  principles  of  simulation 
physics  as  applied  to  sputtering.   This  thesis  is  based 
upon  the  program  QMOL,  a  slightly  modified  version  of  QRAD, 
which  allows  bombardment  by  diatomic,  homonuclear 
molecules.   The  remainder  of  this  section  will  be  presented 
in  terms  of  QRAD.   The  unique  differences  of  QMOL  will  be 
discussed  in  following  sections. 

A  note  on  termionlogy  is  instructive.   Ion  refers  to  a 
single  atom  incident  upon  the  target.   Molecule  refers  to  a 
diatomic,  homonuclear  species  of  the  same  element  as  the 
ion  incident  upon  the  target.   The  molecule  always 
dissociates  upon  impact.   Immediately  following  impact, 
each  fragment  of  the  molecule  is  an  ion.   Atom  refers  to  an 
atom  of  the  monocrystal.   A  sput  is  a  sputtered  atom. 

QRAD  simulates  a  wide  variety  of  experimental 
conditions.   Each  of  the  common  metallic  unit  cells  can  be 
built  and  duplicated  to  form  a  monocrystal  of  up  to  2500 
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atoms  with  a  specific  crystal  plane  exposed  to  the  incident 
ion.   Additionally,  QRAD  can  model  alloys,  stepped 
surfaces,  and  adsorbed  surface  layers. 

The  ion,  of  given  mass,  is  shot  at  the  surface  by 
assigning  to  it  velocity  components  consistent  with  its 
angle  of  incidence  and  energy,  and  positioning  it 
immediately  above  a  specific  impact  point.   Hamilton's 
equations  of  motion  are  solved  numerically  by  an  average 
force  method  [45]  for  the  resultant  MI  collision  cascade, 
or  trajectory. 

The  trajectory  advances  by  timesteps.   During  each 
timestep,  the  positions  and  velocities  of  the  ion  and  atoms 
develop  simultaneously.   Each  ion  and  atom  is  "looked  at" 
and  its  motion  and  forces  analyzed  at  the  end  of  each  time 
step.   The  "look"  each  receives  depends  upon  its  total 
energy  and  the  total  energy  of  its  nearest  neighbors.   A 
hierarchy  of  three  mini-logics  corresponding  to  the  depth 
of  the  "look"  significantly  reduces  computer  time  compared 
with  earlier  generations  of  QRAD.   The  program  flags  all 
ions  and  atoms  which  rise  above  the  surface.   The 
trajectory  terminates  when  the  energy  of  the  most  energetic 
atom  decreases  to  a  specified  value  below  which  the 
probability  of  additional  atoms  rising  above  the  surface  is 
very  small. 

Each  flagged  ion  or  atom  is  tested  by  the  program  to 
determine  whether  it  has  sufficient  outward  momentum  to 
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completely  escape  the  surface  into  the  "vacuum".   Atoms  and 
ions  which  pass  this  test  are  listed  as  sputtered  and 
scattered,  respectively.   The  program  records  the  final 
positions  and  velocities  of  the  sputs  and  scattered  ions 
and  restores  the  system  so  that  another  trajectory  can  be 
executed.   Individual  trajectories  may  be  studied  in  detail 
to  follow  the  sputtering  process  or  identify  the  occurrence 
of  other  phenomena  such  as  channeling. 

Global  results  which  may  be  compared  to  experimental 
data  are  obtained  by  averaging  the  results  over  N 
trajectories,  where  N  is  typically  80-100.   The  sample  of 
trajectories  is  based  upon  N  impact  points  uniformly 
distributed  within  an  impact  area  chosen  such  that  it  is  a 
representative  symmetry  area  of  the  surface.   The  N 
trajectories  determine  sput  yields,  energy  and  angular 
distributions  and  the  spatial  distribution  of  yield  over 
the  impact  area.   Multimer  formation  is  determined  from 
additional  sput  analysis.   Dimers  and  trimers  are 
identified  and  listed  with  their  energy  and  angular 
distributions . 

B.   INPUT  PARAMETERS 

Molecular  data  input  parameters  are  listed  and  defined 
below  so  that  their  meanings  will  be  clear  in  following 
sections . 
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1.   Bullet  parameter  card 
ion  (1) ,  ion  (2 ) 


chemical  symbols  of 
constituent  atoms 


LSS  =  0/3 


flag: 

0  =  single  incident  ion 

3  =  incident  molecule 


THEM 


9m,  polar  angle  of 

ion  (2)  with  respect  to 

ion  (1) 


PHIM 


<j)m,  azimuthal  angle  of 
ion  (2)  with  respect  to 
ion  (1) 


SEP 


p ,  interatomic  separation 
in  angstroms. 


IRANF  =  0/1 


flaq  for  random  molecule 


2. 


ISEED 

Reconstruction  card 
IRECON  =0/1 


seed  value  for  random 
no.  generator 

flag  for  reconstruction 
of  trajectory 


ITRJ 
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no.  of  trajectory  to  be 
reconstructed 


C.   MICROCRYSTAL 

A  360-atom  copper  (FCC)  monocrystal  is  used  for  each 
simulation.   The  crystal  is  "cut"  to  expose  the  (111)  plane 
to  the  incident  molecule  (Fig.  1).   The  crystal  contains  4 
layers  with  90  atoms  per  layer  and  is  designated  (15x4x12): 
15  planes  wide,  4  planes  thick,  and  12  planes  long.   RYBL 
is  the  minimum  positive  distance  an  atom  must  rise  above  a 
surface  before  it  is  considered  ejected.   Note  that  the 
first  numbered  atom  in  the  crystal  is  labelled  atom  3  (Fig. 
2).   Atoms  1  and  2  comprise  the  molecule.   When  comparing 
QMOL  and  QRAD  trajectories,  AT0M(x)QM0L  =  AT0M( x-1 ) QRAD. 

Containment,  the  ability  of  the  target  lattice  to 
contain  the  entire  collision  cascade,  is  discussed 
elsewhere  in  detail  [34].   The  size  of  this  microcrystal  is 
a  trade-off  between  asymptotic  containment  and  reasonable 
computer  run-time.   A  series  of  tests  have  determined  that 
absolute  yield  containment  requires  approximately  150  atoms 
per  layer  at  1.0  kev  and  about  250  atoms  per  layer  at  5.0 
kev.   Absolute  containment  of  an  event  is  paramount  for  an 
analytic  description  of  the  sputtering  yield.   However,  it 
is  less  significant  when  studying  mechanisms  or  observables 
such  as  yield  ratios,  and  energy  and  angular  distributions. 
Experience  shows  [28,-32,  34]  that  this  lattice  provides 
adequate  containment  for  these  observables  at  0.50  to  5.0 
kev  per  ion. 
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D.   ORIENTATION  OF  THE  MOLECULE 

Figure  3  shows  the  orientation  of  the  molecule  with 
respect  to  the  crystal  axes.   Note  that  the  -y  axis  is  the 
outward  normal  for  the  (111)  plane. 

The  molecule  has  5  degrees  of  freedom.  Angles  9  j_  and 
<J>i  define  the  incident  direction  of  the  molecule.  8m  and 
<$>m,  the  polar  and  azimuthal  angles  respectively,  and  the 
interatomic  separations  ,  define  the  position  of  ion(2) 
relative  to  ion(l).  Ion(2)  is  assigned  the  same  velocity 
components  as  ion(l).  Vibrational  and  rotational  degrees 
of  freedom  are  not  included. 

The  ranges  of  the  angles  are  0  _<  9m,  9j_  _<  tt/2 
and   0  j<  <J>m,  <f)i  _<_  2tt.   Therefore,  ion(2)  is  always  trailing 
ion  (1)  prior  to  impact.   The  molecule  can  assume  any 
orientation  with  respect  to  the  surface  because  ion(l)  and 
ion(2)  are  identical.   This  procedure  would  not  be  suitable 
for  an  A3-type  molecule  such  as  CO.   Each  simulation  was 
executed  with  the  molecule  normally  incident  upon  the 
surface.   9m  and  <f>m  were  set  at  specific  values  for  some 
runs,  but  at  random  values  for  others  to  be  more  similar  to 
experimental  conditions.   Some  effects  of  beam  polarization 
can  be  inferred  from  trajectories  in  which  9m  and  <j>m  are 
set  at  specific  values  :   (1)  transversely  polarized 
molecule,  9m  =  tt/2,  0  _<_  cf>m  <_   2tt  ,   (2)  longitudinally 
polarized  molecule,  9m  =  cj>m  =  0. 
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Random  values  of  0  m  and  <f>m  are  generated  as  follows: 
When  IRANF  is  set,  the  program  calls  the  library  random 
number  generator  LRND  which  uses  the  value  of  ISEED  to 
generate  5  random  numbers.   Four  are  used  to  calculate  9 
and  <j)m  by  von  Neumann's  method: 

m 


rn 


(R^  +   r|) 


,      .  -1   2  R,  R, 

<fc   =  sin        3   4 
m 


»3  +    rJ) 


where  0  _<  9  m,  ^m  _<  tt/2.   The  fifth  sets  0  _<  <j>  m  <_  tt  or 

tt  <  <|)  m  j<_  2  tt.   LRND  is  called  for  each  subsequent  trajectory 

and  uses  its  internally  updated  value  of  ISEED  to  generate 

the  next  set  of  random  numbers. 

The  initial  position  of  ion(2)  with  respect  to  the 

impact  point  is  calculated  as  follows: 

DPX    =    p    sin    6m   cos    4>m 

DPY   =   - p   cos  6   m 

DPZ  =  p  sin6  m  sin<J>  m. 

The  exact  impact  point  of  ion(2)  is  unknown  prior  to  impact 
unless  the  molecule  is  transversely  polarized.   In  all 
cases,  ion(2)'s  motion  is  perturbed  by  interaction  with 
ion(l)  at  impact.   Ion(2)'s  impact  point  can  be  determined 
by  examining  the  first  5-10  timesteps  in  the  printed  output 
of  a  trajectory  reconstruction. 

A  random  molecule  trajectory  is  not  reproducible  unless 
the  exact  values  of  8m,  <^m,  and  ISEED  are  repeated  from  the 
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original  trajectory.   A  specific  trajectory  is 
reconstructed  as  follows:   When  IRECON  is  set,  the  program 
calls  LRND  and  executes  a  LRND  loop  ITRJ-times,  where  ITRJ 
is  the  number  of  the  trajectory  to  be  reconstructed.   The 
program  now  holds  the  same  set  of  5  random  numbers  for 
trajectory  ITRJ  as  in  the  original  run  and  calculates  9m 
and  <(>m  as  before. 

E.   IMPACT  AREA 

The  impact  area  for  each  simulation  is  a  rectangle 
containing  a  grid  of  104  uniformly  distributed  impact 
points  (Fig.  4)  and  is  a  representative  symmetry  area  of 
the  entire  (111)  surface  (Fig.  2).   The  impact  area  was 
designed  for  bombarding  a  stepped  Ni  (111)  surface  [46]  and 
was  adapted  to  this  research  to  provide  an  extensive 
sampling  area. 

Note  that  the  reference  target  point  is  centered  on 
atom  36  (for  QMOL)  and  is  designated  (RBX,  RBZ)  =  (6,  4)  in 
units  of  planes  as  measured  from  the  origin  at  the  center 
of  atom  3.   Likewise,  the  impact  points,  (SPX,  SPZ),  are  in 
units  of  planes  and  are  labelled  in  figure  4  with  their 
corresponding  trajectory  numbers.   Figure  5,  drawn  to 
scale,  shows  the  complete  grid  of  impact  points 
superimposed  on  the  surface  with  atom  44  centered  in  the 
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area.   Note  the  positions  of  the  atoms  in  the  second 
layer. 

Independent  data  sets  can  be  generated  from  the 
104-point  data  set  by  simple  displacements  from  the 
original  RBX/RBZ.   Figures  6  and  7  show  two  sets  of  RBX/RBZ 
displacements  which  can  generate  1-12  displaced  data  sets. 
The  utility  of  the  displaced  data  sets  is  to  significantly 
increase  the  sput  population  for  a  particular  simulation 
without  significantly  changing  the  position  of  the  impact 
area.   The  increased  populations  reduce  statistical 
flucuations  in  energy  and  angular  distributions. 

F.   ENERGIES 

Recall  that  Bader,  et .  al.,  [44]  conducted  low-energy 
(<  8  kev/ion)  molecular  sputtering  experiments,  and  that 
investigations  of  non-linear  cascade  effects  have  utilized 
high-energy  (100-200  kev/ion)  molecules.   The  simulations 
presented  in  this  thesis  were  conducted  at  0.50,  1.25, 
2.50,  and  5.00  kev/ion. 

This  is  the  first  investigation  of  molecular  sputtering 
using  the  latest-generation  simulation  model.   An  extensive 
data  base  of  single-ion  simulation  sputtering  experiments 
exists  in  this  energy  range.   A  molecular  sputtering  data 
base  of  similar  scope  will  provide  additional  insight  into 
the  analytic  discrepancies  of  low-energy  sputtering. 
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The  results  of  these  simulations  can  be  compared  with 
the  experimental  results  of  Bader,  et.  al. ,  [44]. 

Andersen  and  Bay  suggest  that  non-linear  yield 
enhancement  is  a  slowly  decreasing  function  of  energy  [39] 
The  trends  of  the  0j/20+  yield  ratios  at  these  low 
energies  will  be  instructive. 

Note  that  the  heat  of  sublimation  of  copper  is  3.53 
ev/atom  and  the  surface  binding  energy  is  2.45 
ev/atom  for  the  potential  functions  set  into  the  model 
[37].   Assume  that  the  incident  ions  (QRAD)  and  molecules 
(QMGL)  deposit  all  their  energy  into  the  crystal.   The 
resulting  energy  per  atom  values  will  bracket  the  heat  of 
sublimation  and  surface  binding  energy  values.   Systematic 
effects  of  the  relative  magnitudes  of  the  energy  density, 
and  the  surface  binding  energy  and  heat  of  sublimation  in 
the  context  of  the  spike  model  [42]  can  be  studied. 

G.   POTENTIAL  FUNCTIONS 

1.   Atom-Atom  potential  functions 

The  program  uses  the  Pot-II  potential  function,  a 
compound  function  which  provides  a  repulsive  Born-Mayer 
(Gibson-2)  "wall  joined  smoothly  to  an  attractive  Morse 
"well"  by  a  cubic  spline  at  the  points  RA  =  0.83  LU  and 
RB  =  1.10  LU.   See  figure  8. 
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The  repulsive  wall  dominates  the  collision  cascade 
dynamics.   Attractive  forces  have  little  effect  upon  the 
dynamics  but  dominate  the  atom's  ability  to  escape  the 
surface. 

The  potential  function  well  depth  reproduces  the 
solid's  heat  of  sublimation  (total  energy)  by  summing  ever 
all  possible  pairwise  interactions.   The  Cu-Cu  potential 
includes  only  nearest  neighbor  and  next-nearest  neighbor 
interactions  because  it  is  truncated  at  a  radius  Rc  = 
2.40  LU.   The  well  also  reproduces  the  solid's  static 
surface  binding  energy  (SBE)  at  its  pre-tra jectory  level. 
The  SBE  is  the  attraction  a  sput  must  overcome.   The  values 
of  the  heat  of  sublimation  and  SBE  are  3.53  ev/atom  and 
2.45  ev/atom,  respectively. 

See  reference  [38]  for  a  detailed  discussion  of  the 
atom-atom  potential  function. 

2.  Ion-Atom  Potential  Function 

The  Cu/0+  ion-atom  potential  function  is  a 
compound  functon  consisting  of  a  0.7-Moliere  wall  splined 
to  a  Morse  well  at  RA  =  0.44  LU  and  RB  =  0.55  LU, 

Refer  to  reference  [37]  for  a  detailed  discussion 
of  the  ion-atom  potential  function. 

3.  Ion-Ion  Potential  Function 

The  oxygen  molecule  potential  function  is  a  Morse 
potential  truncated  at  Rc  =  2.40  LU  (Fig.  8). 
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IV.   DISCUSSION  AND  RESULTS 

A.   ENERGY  DISTRIBUTIONS 

The  characteristics  of  the  sputtered  atom  energy- 
distributions  for  the  incident  energies  of  0.5,  1.25,  2.5, 
and  5.0  kev/ion  are  compared  for  the  following  simulations: 
(1)  QRAD,  (2)  random  molecule,  QMOL(R),  and  (3)  longitu- 
dinally polarized  molecule,  QMOL(LP) .   These  distributions 
are  equivalent  to  those  which  would  be  obtained  from  an 
ideal  experiment  which  bombards  a  clean,  monocrystalline 
(111)  copper  surface  in  a  perfect  vacuum,  and  collects  and 
analyzes  each  sputtered  atom. 

1 .   Normalized  Energy  Distributions,  0-20  ev/atom 

Figures  9-12  are  normalized  energy  distribution 
curves  comparing  QRAD,  QMOL(R),  and  QMOL(LP)  for  each 
incident  energy  and  a  sputtered  atom  energy  range  of  0-20 
ev/atom.   These  curves  were  obtained  from  total  sputtered 
atom  populations  of  700-1800  atoms  and  are  plotted  as  raw 
data.   Individual  sputtered  atom  populations  are  indicated 
with  each  distribution.   These  small  populations  produce 
distributions  with  many  structural  details  which  are  not 
statistically  significant.   The  smooth  curves  fitted  to  the 
data  have  no  physical  significance  for  any  of  the 
distributions  taken  alone.   However,  when  viewed  in  a 
series,  certain  structural  trends  may  be  elucidated  for 
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specific  simulations  at  different  incident  energies  or 
different  simulations  at  a  specific  energy.   Several 
features  of  the  distributions  are  strikingly  apparent.   The 
maxima  occur  at  precisely  the  same  energy,  approximately 
1.8  ev/atom.   The  only  exception  is  the  maximum  for  QRAD  at 
0.5  kev/ion.   However,  its  maximum  at  approximately  2.2 
ev/atom  is  within  the  statistical  uncertainty  of  the 
simulation  and  the  difference  represents  a  flucuation  of 
8-9  sputtered  atoms. 

A  shoulder  appears  in  each  distribution  between  2 
and  5  ev/atom.   The  shoulders  for  QRAD  and  QMOL(LP) 
decrease  in  energy  and  increase  in  fraction  of  total  yield 
as  incident  energy  increases.   The  shoulders  for  QMOL(R) 
remain  nearly  constant  as  incident  energy  increases:   The 
widths  of  the  shoulders  are  1-1.5  ev/atom  and  the  yield 
fractions  are  0.040-0.054.   This  is  equivalent  to  a  yield 
flucuation  of  41-95  sputtered  atoms.   Secondary  shoulders 
may  be  present  in  the  6-10  ev/atom  range  as  well  as  valid 
structure  at  higher  energies,  but  the  populations  are 
insufficient  to  define  additional  features. 

2.   Normalized  Energy  Distributions,  0-3  ev/atom 

It  is  instructive  to  study  the  normalized  energy 
distributions  for  the  sputtered  atom  energy  range 
containing  the  peaks  seen  in  the  preceding  section. 
Figures  13-16  are  distributions  for  the  energy  range  0-3 
ev/atom.   These  distributions  are  plotted  from  a  "sliding-4" 
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average  of  the  raw  data  as  follows:   raw  data  points  1-4 
are  summed,  followed  by  points  2-5,  3-6,  4-7,  etc.   The  set 
of  averaged  points  is  normalized  and  plotted  as  the 
fraction  of  total  yield  vs.  energy  per  sputtered  atom.   As 
before,  the  smooth  curves  fitted  to  the  data  have  no 
individual  physical  significance. 

The  relative  locations  of  the  apparent  maxima  of 
the  distributions  are  consistently  skewed  toward  lower 
energies  as  indicated  by  lines  A  in  the  figures. 
Significant  shifts  of  the  distributions  for  a  particular 
simulaton  as  a  function  of  incident  energy  are  not  seen. 
In  fact,  the  constant  energy  peaks  of  figures  9-12  suggest 
that  no  significant  shifts  should  occur  in  the  0-3  ev/atom 
range . 

Several  apparent  peaks  are  seen  in  the  distribu- 
tions.  Consider  the  QRAD  distributions.   The  peaks  are 
fitted  credibly  because  of  their  consistency  as  the 
incident  ion  energy  increases.   The  peak  vanishes  at  5.0 
kev/ion,  but  note  that  the  total  sputtered  atom  population 
at  5.0  kev/ion  is  only  half  of  the  populations  at  the  lower 
incident  energies.   The  5.0  kev/ion  QRAD  simulation 
consists  of  1352  trajectories.   Nearly  3000  trajectories 
are  necessary  to  produce  a  sputtered  atom  population 
sufficiently  large  to  verify  a  peak.   No  comparable  peak 
behavior  is  seen  for  QMOL(R)  or  QMOL(LP). 
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Several  low-energy  shoulders  may  be  real  and 
verifiable  for  significantly  increased  populations. 

Assume  that  these  structural  features  are  real. 
Individual  sets  of  trajectories  isolating  the  atoms 
comprising  the  features  can  be  studied  to:   (1)  identify 
the  sputtering  mechanism  generating  the  QRAD  peaks  and 
shoulders,  and  (2)  elucidate  the  molecular  effects 
responsible  for  breaking  up  or  reinforcing  the  mechanisms. 

3 .   Incident  Energy  Dependence  of  the  Normalized  Energy 
Distributions 

Figures  17-22  show  the  ion/molecule  energy 
dependence  of  the  normalized  energy  distribution  functions 
for  each  simulation.   The  distributions  are  identical  to 
those  previously  studied  and  are  displaced  vertically  along 
an  incident  energy  scale.   Fraction  of  total  yield  cannot 
be  read  from  this  scale  but  an  increment  of  fractional 
yield  is  indicated  for  each  figure. 

The  0-20  ev/atom  distributions,  figures  17-19, 
clearly  show  the  alignment  of  the  peaks  and  the  shifts  in 
shoulder  positions.   The  QRAD  shoulders  shift  linearly  to 
lower  energy  and  higher  fraction  as  incident  energy 
increases.   The  QMOL(R)  shoulders  remain  constant  in  energy 
with  no  significant  shift  in  yield  fraction.   The  QMOL(LP) 
shoulders  show  a  non-linear  shift  to  lower  energy  and 
higher  fraction. 
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Figures  20-22  are  the  0-3  ev/atom  distributions. 
The  QRAD  peaks  remain  in  the  vicinity  of  2  ev/atom.   The 
prospective  location  of  a  peak  in  the  5.0  kev/ion 
distribution,  discussed  in  the  preceding  section,  is  clear 
in  figure  20.   No  significant  trends  in  fine  structure  are 
evident  in  the  QMOL(R)  or  QMOL(LP)  distributions.   Real 
shifts  in  the  distributions  for  either  simulation,  as  a 
function  of  incident  energy,  are  not  supported. 
4.   Cumulative  Energy  Distribution  Functions 

Figures  23-25  are  the  cumulative  energy  distribu- 
tion functions  (CDF)  which  plot  the  fraction  of  sputtered 
atoms  with  energy  greater  than  the  indicated  energy  as  a 
function  of  sputtered  atom  energy  for  the  range  0-200 
ev/atom.   The  figures  compare  the  simulations  as  follows: 
(1)  Fig.  23,  QRAD  vs.  QMOL(R),  (2)  Fig.  24,  QRAD  vs. 
QMOL(LP),  and  (3)  Fig.  25,  QMOL(R)  vs.  QMOL(LP).   The 
incident  ion  energy  is  indicated  with  each  curve.   Note 
that  for  a  specific  incident  ion  energy  the  lower  of  the 
two  curves  shows  a  shift  to  the  low  energy  region  of  the 
range  since  each  plotted  point  represents  the  fraction 
greater  than  the  energy  indicated. 

Figure  2  5  shows  the  only  consistent  trend  for  the 
distributions.   The  QMOL(R)  sputtered  atom  populations  are 
shifted  to  lower  energies  than  those  for  QMOL(LP). 
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Figure  23  will  be  particularly  instructive  in 
correlating  the  angular  distribution  effects  discussed  in 
the  following  section. 

B.   ANGULAR  DISTRIBUTIONS 

Figures  26-29  are  the  angular  distributions  of 
sputtered  atoms  (spot  patterns)  for  QRAD  and  QMOL(R). 
These  figures  are  equivalent  to  spot  patterns  obtained 
experimentally  by  collecting  sputtered  deposits  on  flat 
glass  plates. 

The  most  obvious  feature  of  the  spot  patterns  is  the 
six-fold  symmetry  which  has  been  shown  experimentally  [14] 
to  be  a  characteristic  of  preferential  ejection  from  the 
(111)  surface.   The  diffuse  distribution  between  the 
six-fold  axes  is  caused  by  low-energy  (<20  ev)  sputtered 
atoms  ejected  late  in  the  trajectories  when  extensive 
surface  disruption  has  broken  down  the  preferred  ejection 
mechanisms  [47]. 

The  density  of  the  spot  patterns  is  directly  related  to 
the  sputtered  atom  population.   However,  there  are  no 
consistent  trends  for  QRAD  vs.  QMOL(R)  as  a  function  of 
incident  energy.   The  QMOL(R)  spot  patterns  are  not  always 
"darker"  or  more  diffuse  than  the  corresponding  QRAD  spot 
patterns . 

The  relative  densities  of  the  spot  patterns  for  QRAD 
and  QMOL(R)  for  each  incident  energy  can  be  correlated  to 
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figure  23,  the  CDF  for  QRAD  vs.  QMOL(R).   Consider  the 
0-20  ev/atom  energy  range  for  each  incident  energy. 

1.  0.5  kev/ion 

The  six-fold  axes  for  QRAD  in  figure  26a  are  very 
well  defined  but  the  spot  pattern  for  QMOL(R),  figure  26b, 
is  heavily  filled  between  the  axes.   Figure  23a  shows  that 
the  CDF  for  QMOL(R)  is  shifted  more  toward  lower  energies 
than  for  QRAD.   The  higher  fraction  of  low-energy  sputtered 
atoms  is  responsible  for  the  heavy,  diffuse  distribution 
between  the  axes . 

2.  1.25  kev/ion 

The  CDF  for  1.25  kev/ion  incident  energy,  figure 
23b,  is  exactly  opposite  that  for  0.5  kev/ion  but  the 
differences  in  fraction  are  not  as  great.   This  is 
reflected  in  the  spot  patterns,  figures  27a  and  27b.   Each 
is  very  dark  but  the  QRAD  pattern  now  shows  slightly  more 
filling  between  the  axes  than  the  QMOL(R)  pattern. 

3 .  2 . 50  kev/ion 

The  CDF  and  spot  patterns  for  2.50  kev/ion  incident 
energy,  figures  23c  and  28,  respectively,  follow  the  same 
trend  as  those  for  1.2  5  kev/ion  incident  energy.   The 
fractional  differences  in  the  CDR  are  smaller  but  the 
higher  QRAD  sputtered  atom  population  produces  a  darker 
spot  pattern  (  the  percentage  difference  in  sputtered  atom 
population  between  QRAD  and  QMOL(R)   is  greatest  at  this 
incident  energy) . 

35 


4.   5.0  kev/ion 

The  CDF  curves  for  5.0  kev/ion  incident  energy, 
figure  23d,  are  virtually  identical,  and  the  spot  patterns 
for  this  energy,  figures  29a  and  29b,  are  nearly 
indistinguishable. 

C.   SPUTTERING  YIELD  RATIOS 

Figure  30  compares  the  experimental  yield  ratios,  for 
the  sputtering  of  polycrystalline  copper  by  N+  and 
N+2  [44],  and  the  simulation  sputtering  yield  ratios. 
The  sputtering  yield  ratio  is  defined  as  follows: 
ratio  =  yield  from  X2+ 


2 (yield  from  X+) 

where  X  is  the  bombarding  element.   Two  sets  of  simulations, 
"real"  and  "artificial",  are  considered. 

1.  Experimental 

The  experimental  sputtering  yield  ratios  for  N+ 
and  N2+  incident  upon  polycrystalline  copper  [44]  are 
plotted  in  figure  30a  and  30b  for  the  incident  energy  range 
0-8  kev/ion.   A  very  good  straight-line  fit  is  obtained 
which  shows  a  slowly  increasing  yield  ratio  with  increasing 
incident  energy. 

2 .  "Real"  Simulations 

The  sputtering  yield  ratios  for  the  "real" 
simulations,  QMOL(R)  and  QMOL(LP),  are  plotted  in  figure 
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30a.   The  connecting  lines  have  no  physical  significance. 
In  general,  the  simulation  yield  ratios  increase  much  more 
slowly  with  incident  energy  than  the  experimental  values. 
The  ratios  for  QMOL(LP)  are  consistently  lower  than  those 
for  QMOL(R)  and  both  dip  considerably  at  1.25  kev/ion 
incident  energy. 

3.   "Artificial"  Simulations 

The  sputtering  yield  ratios  for  two  "artificial" 
simulations  are  plotted  in  figure  30b.   In  the  first  case, 
Ar+  and  Ar^  are  substituted  for  0+  and  O2"1" 
but  the  Cu-0  ion-atom  potential  is  retained.   Therefore, 
the  mass  ratios  are  changed  for  QRAD  and  QMOL(R)  from  0.251 
and  0.503  for  0+/02+  to  0.628  and  1.257  for 
Ar+/Ar2+/  respectively.   In  the  second  case,  the 
B-potential  [37]  is  substituted  for  the  Cu-0  potential  and 
the  oxygen  mass  is  retained.   The  ratios  for  ArJ/Ar* 
are  lower  than  the  corresponding  "real"  QMOL(R)  ratios  in 
figure  30a  and  an  identical  dip  occurs  at  1.25  kev/ion  is 
almost  totally  absent. 

D.   EFFECTS  OF  MOLECULAR  ORIENTATION 

The  effects  of  molecular  orientation  upon  the 
sputtering  yield  ratio  are  examined  for  the  following:   (1) 
longitudinally  polarized  molecule,  QMOL(LP),  (2) 
transversely  polarized  molecule  (9 m  =  90  deg.),  QMOL(TP) 
and  (3)  polar  angle  8m  =  45  deg.,  QMOL(45).   Figures  31a 
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and  31b  show  the  sputtering  yield  ratio  as  a  function  of 
azimuthal  angle,  <j>m,  at  1.25  kev/ion  incident  energy  for 
QMOL(TP)  and  QMOL(45).   Lines  of  constant  yield  ratio  are 
indicated  for  QMOL(LP)  at  1.25  kev/ion  incident  energy  and 
the  average  QMOL(LP)  yield  ratio  for  the  four  incident 
energies.   Note  that  azimuthal  angle  is  not  defined  for  a 
longitudinally  polarized  molecule. 
1 .   General 

On  the  average  for  a  random  molecular  orientation, 
the  two  dissociated  ions  will  not  collide  with  the  surface 
simultaneously.   When  the  molecule  is  transversely  polar- 
ized, the  time  separation  between  collisions  is  minimized 
and  the  separation  between  impact  points  is  the  greatest. 
The  degree  of  overlap  of  the  individual  collision  cascades 
will  be  determined  by  the  separation  in  the  diatomic 

o 

molecule,  1.22  A  for  02-   When  the  molecule  is 
longitudinally  polarized,  the  time  separation  is  maximized 
and  the  impact  points  are  very  nearly  identical. 

The  longitudinally  polarized  molecule  produces  the 
highest  energy  density  in  the  collision  cascade  since  the 
energy  from  both  ions  is  deposited  in  the  same  region  of 
the  lattice.   If  energy  density  were  the  controlling 
factor,  the  longitudinally  polarized  molecule  should 
produce  the  highest  yield  ratios. 

The  trailing  ion  in  the  longitudinally  polarized 
molecule  is  colliding  with  a  surface  already  disrupted  by 
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the  leading  ion  and,  in  general,  will  not  "see"  a  surface 
layer  atom  before  impact.   In  adition,  the  trailing  ion 
overtakes  the  leading  ion  and  is  deflected  slightly  by  the 
leading  ion's  collision  with  the  surface.   These  factors 
may  result  in  increased  penetration  depths  and  a  higher 
incidence  of  channeling  for  bombardment  by  longitudinally 
polarized  molecules. 

2.  Yield  Ratios 

The  yield  ratios  for  QMOL(LP)  and  QMOL(45)  at  1.25 
kev/ion  incident  energy  are  plotted  in  figure  31,  with 
QMOL(LP)  yield  ratios  indicated,  as  a  function  of  the 
azimuthal  angle  4>m.   The  ratios  for  both  cases  show  a  good 
linear  fit  and  are  nearly  identical.   The  QMOL(LP)  values 
are  not  significantly  different  from  either  QMOL(TP)  or 
QMOL(45) . 

3 .  Scattered  ions 

The  following  table  summarizes  and  compares  total 
yields  and  the  number  of  scattered  ions  for  QRAD,  QMOL(R) 
and  QMOL(LP) : 

no. 


tra  j 

total  yield 

scattered  ions 

0. 5  kev/ion 

QRAD 

QMOL(R) 

QMOL(LP) 

936 
520 
520 

1665 
1770 
1657 

104 
225 
138 

1.25  kev/ion 

QRAD 

QMOL(R) 

QMOL(LP) 

936 
520 
520 

1653 
1660 
1437 

72 
67 
42 

2 . 50  kev/ion 

QRAD 

QMOL(R) 

QMOL(LP) 

1352 
520 
520 

1439 
1181 
1069 

20 
29 

15 

5.0  kev/ion 

QRAD 

QMOL(R) 

QMOL(LP) 

1352 
520 
520 

927 
760 
746 

8 
6 
2 
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E.   OTHER  MOLECULAR  EFFECTS 
1 .   Yield  per  Layer 

Detailed  analyses  of  single-ion  sputtering 
trajectories  have  shown  that  the  majority  of  sputtered 
atoms  come  from  sputtering  mechanisms  which  are  completely 
confined  to  the  surface  and  near-surface  layers  [29]. 
Consequently,  the  majority  of  sputtered  atoms  come  from  the 
surface  layer. 

The  sputtering  yields  for  each  layer  of  the 
microcrystal  are  listed  below,  for  each  simulation  and 
incident  energy,  with  the  fractional  yields  from  layer  2 
(L2/N)  and  all  subsurface  layers  (  (L2  +  L3  +  L4.)/N) 
indicated. 


Total  yield, 

#  Each  Layer 

hz 

L2+L3+L4 

0 . 5  kev/ion 

N 

L1/L2/L3/L4 

N 

N 

QRAD 

1665 

1635/28/2/0 

0.016 

0.018 

QMOL(R) 

1770 

1708/53/7/2 

0.029 

0.035 

QMOL(LP) 

1657 

1589/54/13/1 

0.032 

0.041 

1.25  kev/ion 

QRAD 

1653 

1626/23/2/2 

0.013 

0.016 

QMOL(R) 

1660 

1614/31/15/0 

0.018 

0.027 

QMOL(LP) 

1437 

1397/35/5/0 

0.024 

0.027 

2 . 50  kev/ion 

QRAD 

1439 

1420/16/3/0 

0.011 

0.013 

QMOL(R) 

1181 

1147/31/3/0 

0.026 

0.028 

QMOL(LP) 

1069 

1033/35/1/0 

0.032 

0.033 

5.0  kev/ion 

QRAD 

927 

909/17/1/0 

0.018 

0.019 

QMOL(R) 

760 

737/21/2/0 

0.027 

0.030 

QMOL(LP) 

746 

718/26/2/0 

0.034 

0.037 
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Note  that  for  molecular  sputtering,  the  vast  majority  of 
sputtered  atoms  still  originate  in  the  surface  layer. 
However,  the  relative  yields  from  subsurface  layers  are 
much  greater  for  molecular  than  for  single-ion  sputtering. 
2.   Normalized  Sputtering  Frequency 

Figures  32-39  are  the  surface  layer  sputtering 
frequency-location  diagrams  for  QRAD  and  QMOL(R)  at  each 
incident  energy.   Each  diagram  is  normalized  to  the  most 
frequently  sputtered  atom  for  the  simulation:   a  circle 
containing  ten  rings  indicates  the  most  frequently 
sputtered  atoms ( s ) .   A  circle  containing  zero  rings  is  an 
atom  that  never  sputtered  and  one  ring  indicates  an  atom 
atom  that  sputtered  at  least  once  in  the  entire  set  of 
trajectories . 

Atom  43  (44),  for  QRAD  (QMOL),  is  the  primary 
knock-on  atom  (PKA)  in  at  least  46%  of  the  trajectories 
(see  figure  5).   It  is  always  one  of  the  least  sputtered 
atoms.   In  fact,  its  frequency  decreases  as  incident  energy 
increases.   Note  that  the  PKA  sputtering  frequencies  are 
always  equal  for  QRAD  and  QMOL(R)  at  each  incident  energy, 
with  the  exception  of  2.50  kev/ion:   at  this  energy  the 
frequencies  differ  by  less  than  0.1. 

At  0.5  kev/ion,  the  most  frequently  sputtered  atoms 
are  both  nearest  and  next-nearest  neighbors,  relative  to 
the  PKA.   But  at  1.25,  2.5,  and  5.0  kev/ion  the  most 
frequently  sputtered  atoms  are  almost  exclusively  next  to 
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the  most  frequently  sputtered  atoms  are  triangular  or 
diamond-shaped. 

F.   THE  MODEL 

Computer  run-times  for  QMOL  are,  on  the  average,  nearly 
twice  those  for  QRAD  at  the  same  energy  per  incident  ion. 
This  is  reasonable  since  a  QMOL  simulation  is  essentially  a 
superposition  of  two  single-ion  collision  cascades.   Run- 
time increases  as  incident  energy  decreases  because  of  ion 
size  effects.   Several  QMOL  simulations,  particularly  at 
0.50  kev/ion,  had  to  be  split  into  two  separate  runs  to 
avoid  exceeding  the  computer  system's  job  time  constraints. 

The  reproducibility  of  a  QMOL  trajectory  is  extremely 
sensitive  to  the  precision  of  9m  and  <J>m.   The  IBM  360 
system  is  designed  to  retain  7  significant  digits,  in 
single-precision,  and  this  is  the  precision  of  the  random 
angles  9m  and  <j>m.   An  attempt  to  reconstruct  a  trajectory 
by  manually  entering  9m  and  $rn  to  3  significant  digits 
produced  only  2  of  the  original  10  sputtered  atoms  for  the 
trajectory.   Subsequent  entries  of  4  and  5  significant 
digits  produced  yields  of  4  and  7  sputtered  atoms, 
respectively.   The  reconstruction  loop  reproduces  a 
trajectory  exactly. 
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V.   CONCLUSIONS 

A.  THE  MODEL 

This  research  reinforces  the  feasibility  and 
flexibility  of  computer  simulated  sputtering.   The  single- 
ion  sputtering  model  was  easily  modified  to  simulate 
homonuclear,  diatomic  molecular  sputtering  within  the 
original  constraints  and  basic  physical  assumptions.   No 
additional  approximations  or  rethinking  of  the  physics 
inherent  in  the  basic  model  were  necessary. 

B.  SPUTTERED  ATOM  NORMALIZED  ENERGY  DISTRIBUTIONS 
Recall  the  alignment  of  the  peaks  in  the  0-20 

ev/sputtered  atom  distributions.   The  peaks  are  clearly  not 
an  incident  energy,  energy  density,  or  molecular  orienta- 
tion effect.   Preliminary  results  of  other  simulations,  not 
included  here,  show  that  the  peaks  are  insensitive  to  the 
mass  ratio,  but  do  shift  to  lower  energies  when  the 
ion-atom  potential  function  is  changed  from  the  Cu-0+ 
compound  potential  to  the  B-potential.   The  shift  is  nearly 
identical  for  single-ion  and  molecular  sputtering. 

The  0-3  ev/sputtered  atom  energy  distributions  provide 
the  greatest  insight  into  the  effects  of  the  molecule  upon 
the  sputtering  process.   The  QRAD  peaks  for  0.50,  1.25,  and 
2.50  kev/ion  incident  energy  may  be  real.   A  mechanism 
unique  to  the  QMOL  cascades  is  clearly  breaking  up  the 
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peaks  and  in  one  case  (QMOL(R)  at  2.50  kev/ion)  may  be 
reinforcing  the  peak.   The  source  of  this  molecular 
sputtering  mechanism  may  be  linked  to  the  substantially 
increased  fraction  of  sputtered  atoms  from  subsurface 
layers  for  the  QMOL  simulations. 

C.   YIELD  RATIOS 

The  simulation  sputtering  yield  ratios  agree  reasonably 
well  with  experimental  data.   Two  primary  factors  may  cause 
the  differences  in  yield  ratios:  (1)  the  difference  in  ion- 
atom  potential  functions,  and  (2)  the  polycrystalline 
target  used  in  the  laboratory  experiments.   Mass  ratio  is 
probably  not  significant  since  the  mass  ratios  for  N  and  0 
differ  by  only  12%. 

The  majority  of  the  yield  ratios  are  less  than  one  and, 
in  general,  the  overall  trend  for  the  yield  ratios  is  a 
slow  increase  with  increasing  incident  energy.   Energy 
density  is  not  an  accurate  predictor  of  the  sputtering 
yield  ratio.   In  fact,  the  yield  ratios  for  the 
longitudinally  polarized  molecule,  the  configuration 
producing  the  highest  energy  density  in  the  lattice,  are 
all  lower  than  the  ratios  for  the  corresponding  random 
molecule  simulations. 

The  anomalous  dip  in  the  sputtering  yield  ratios  at 
1.25  kev/ion  may  be  an  artifact  of  the  ion-atom  potential 


44 


function  since  it  vanishes  in  the  B-potential  "artificial" 
simulation.   It  is  interesting  to  note  that  1.25  kev 
averaged  over  the  360  -  atom  microcrystal  is  only  0.06 
ev/atom  less  than  the  heat  of  sublimation  of  the 
microcrystal  modelled  by  the  Cu-Cu  potential  function. 

D.   MOLECULAR  ORIENTATION 

Sputtering  yield  ratio  is  strikingly  insensitive  to 
molecular  orientation  at  1.25  kev/ion  incident  energy. 
Although  QMOL(TP)  and  QMOL(45)  yield  ratios  are  not 
available  for  the  other  incident  energies,  no  significant 
difference  is  expected  since  the  QMOL(R)  and  QMOL(LP) 
ratios  are  close  at  those  energies. 

Molecular  beam  polarizability  may  have  a  definite 
application  for  ion  implantation.   The  QMOL(LP)  simulations 
produced  substantially  less  scattered  incident  ions  at  • 
1.25,  2.5,  and  5.0  kev/ion. 
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VI.   RECOMMENDAT IONS 

A.  THE  MODEL 

QMOL  can  easily  be  modified  to  simulate  sputtering  by 
an  AB-type  molecule  as  follows.   Change  the  random  number 
generator  loop  to  produce  six  random  numbers  instead  of 
five.   Use  the  sixth  random  number  to  set  0  <  8  m  _<_  tt  or  tt  < 
9m  _<_  2tt  in  the  same  way  that  the  fifth  random  number  sets 
the  value  of  $m.       The  atom  that  is  in  the  lower  hemisphere 
of  the  molecule's  locus  of  orientation  becomes  ion(l).   The 
program  then  sets  up  the  molecule  at  the  impact  point.   In 
principle,  the  program  can  be  modified  to  simulate 
sputtering  by  triatomic  or  higher-order  polyatomic 
molecules . 

B.  SPUTTERED  ATOM  NORMALIZED  ENERGY  DISTRIBUTIONS 

It  is  necessary  to  identify  specific  atoms  comprising 
the  QRAD  peaks.   Then,  individual  QMOL  trajectories  common 
to  these  atoms  can  be  examined  in  detail  to  elucidate  how  a 
peak  is  dispersed  into  a  much  wider  energy  band.   However, 
sputtered  atom  populations  for  QRAD  at  5.0  kev/ion  and  all 
QMOL  simulations  are  too  sparse  and  statistically  uncertain 
for  this  detailed  analysis.   The  5.0  kev/ion  QRAD 
simulation  should  be  run  with  a  new  displaced  data  set 
designed  to  augment  the  1352  trajectories  already  executed 
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and  produce  a  sputtered  atom  population  sufficiently  large 
to  confirm  or  negate  the  appearance  of  a  peak.   Likewise, 
each  QMOL  sputtered  atom  population  should  be  augmented  by- 
runs  with  the  displaced  data  set  no.  2  impact  points.   The 
combination  of  much  larger  sputtered  atom  populations  and 
the  sliding-4  averaging  technique  will  enable  more  precise 
correlation  between  structural  features  in  the  0-3 
ev/sputtered  atom  energy  distributions. 

C.   YIELD  RATIOS 

A  systematic  study  of  the  differences  between  the 
simulation  and  experimental  sputtering  yield  ratios  must 
focus  upon  variations  in  the  ion-atom  potential  function 
until  laboratory  data  for  single-crystal  molecular 
sputtering  becomes  available.   Two  primary  options  are 
available:   (1)  Model  the  Cu-N+  potential  while  retaining 
the  oxygen  mass  in  the  simulation,  and  (2)  increase  the 
O2  bond  length  to  reduce  cascade  overlap. 

Individual  trajectories  should  be  compared  for  QMOL(R) 
at  0.50  and  1.25  kev/ion,  and  QMOL(LP)  at  1.25  and  2.50 
kev/ion  to  bracket  the  dip  at  its  largest  variation.   The 
relative  positions  of  the  ions  for  each  type  of  simulation 
can  be  traced  from  timestep  to  timestep  in  the  individual 
trajectories.   This  may  show  peculiar  differences  in 
cascade  overlap  causing  inconsistent  yield  ratios. 
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D.   ION-ATOM  POTENTIAL  FUNCTION 

The  cubic  spline  in  the  Cu-0+  potential  function  must 
be  adjusted.   A  slight  knee  in  the  curve  occurs  at  approx- 
imately 100  ev  potential  energy  and  0.5  LU  separation  when 
the  function  is  plotted  from  the  program's  potential  table. 
This  shows  that  the  spline  is  not  precisely  matching  the 
magnitude  and  slope  of  the  Moliere  and  Morse  potentials  at 
the  current  values  of  RA  =  0.44  LU  and  RB  =  0.55  LU . 
For  an  initial  iteration  change  RA  and  Rg  to  0.40  and 
0.60  LU,  respectively,  execute  the  program  for  one  timestep 
to  print  a  new  potential  table,  and  replot  the  ion-atom 
potential  function.   A  change  in  the  spline  will  not  affect 
the  simulation  results. 
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Fig.  3.   Molecule  orientation  with  respect 
to  the  (111)  surface  of  the  target 
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Fig.  5.  Impact  area  on  the  (111)  surface,  drawn  to 
scale.  The  central  atom  is  atom  44  (QMOL) 
The  dashed  circles  are  second-laver  atoms. 
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Fig.  23.   Cumulative  energy  distributions  of  sputtered  atoms. 
C  +  QRAD,  o  QMOL(R)) 
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Fig.  24.   Cumulative  energy  distributions  of  sputtered  atoms 
(  +  QRAD,  o  QMOL(LP)) 
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Fig.  25.   Cumulative  energy  distributions  of  sputtered  atoms 
(  +  QMOL(R),  o  QMOL(LP)) 
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Fig.  26.   Spot  patterns,  0.50  kev/ion, 
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Fig.  27.   Spot  patterns,  1.25  kev/ion. 
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Fig.  28.   Spot  patterns,  2.50  kev/ion, 
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Fig.  29.   Spot  patterns,  5.00  kev/ion. 
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Fig.  30.   Experimental  [44]  and  simulation  sputtering  yield 
ratios  vs.  incident  energy. 
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Fig.  31.   Simulation  sputtering  yield  ratio  vs.  azimuthal  angle,  <f>  , 
for  QMOL(TP)  and  QMOL(45). 

(Line  (a)  is  the  QMOL(LP)  avg.  yield  ratio;   line  (b)  is 
the  QMOL(LP)  yield  ratio  at  1.25  kev/ion.) 
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Fig.  32.   QRAD  normalized  sputtering  frequency,  0.50  kev/ion, 
See  text  for  explanation  of  scaling. 
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Fig.  33.   QMOL(R)  normalized  sputtering  frequency,  0.50  kev/ion, 
See  text  for  explanation  of  scaling. 
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Fig.  34.   QRAD  normalized  sputtering  frequency,  1.25  kev/ion, 
See  text  for  explanation  of  scaling. 
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Fig.  35.  QMOL(R)  normalized  sputtering  frequency,  1.25  kev/ion, 
See  text  for  explanation  of  scaling. 
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Fig.  36.   QRAD  normalized  sputtering  frequency,  2.50  kev/ion, 
See  text  for  explanation  of  scaling. 
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Fig.  37.   QMOL(R)  normalized  sputtering  frequency,  2.50  kev/ion, 
See  text  for  explanation  of  scaling. 
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Fig.  38.   QRAD  normalized  sputtering  frequency,  5.00  kev/ion. 
See  text  for  explanation  of  scaling. 
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Fig.  39.   QMOL(R)  normalized  sputtering  frequency,  5.00  kev/ion, 
See  text  for  explanation  of  scaling. 
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